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Abstract. The dynamics of the surface plasmon in laser-driven metallic nanoparticles is described by means 
of a master-equation formalism. Within the Markov approximation, the dynamics is studied for different 
regimes ranging from weak excitation in photoabsorption experiments to strong excitation in pump-probe 
spectroscopy. It is shown that two collective levels are sufficient to describe the dynamics of the surface 
plasmon. On this basis, we predict the appearance of sidebands in the absorption spectrum of the probe 
laser field in pump-probe experiments. 
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1 Introduction 

A surface plasmon excited by a pump laser in a metallic 
nanoparticle decays on the scale of tens of femtoseconds 
by heating the electron system [1,2]. The equilibration 
with the ionic lattice occurs on the much longer timescale 
of a few picoseconds. This electronic relaxation has been 
thoroughly studied by means of pump-probe experiments 
[3-6]. Present femtosecond resolution experiments allow 
to concentrate on the initial dynamics of the electronic 
degrees of freedom [7-9]. While the dynamics on the pi- 
cosecond scale has been analysed with quasi-equilibrium 
theories assuming the thermalisation of the electron gas as 
a whole [10], much less is known for the initial dynamics. 
It is precisely this initial regime on which we focus in this 
work. 

For the description of the initial dynamics of the elec- 
tron gas, before it has thermalised after an excitation by 
a laser pulse, it is appropriate to decompose the electron 
system into two parts [11-13]. The main degree of freedom 
describes the surface plasmon and refers to the centre of 
mass of the electron gas. The remaining degrees of free- 
dom consist of the relative coordinates. If their number 
is sufficiently large, this electronic subsystem behaves to 
a good approximation as a heat bath to which one can 
attribute an effective temperature. 

The surface plasmon coupled to the relative coordi- 
nates and driven by an external field is reminiscent of the 
situation encountered in quantum optics, where a laser- 
driven atom is coupled to an environment consisting of 
the electromagnetic field modes [14]. This analogy leads 
us to predict the appearance of sidebands in the light ab- 



sorption of a weak probe beam while the nanoparticle is 
illuminated by a strong pump laser. This effect should be 
experimentally observable and offers a possibility to study 
the surface plasmon dressed by photons of the pump laser 
field. 

We start in the next section by introducing the model 
of a laser-driven nanoparticle. In Section [3l a master- 
equation description for the dynamics of the electronic 
centre of mass coupled to the relative-coordinate system 
is developed. We will argue that the Markov approxima- 
tion applies in our situation. In order to make progress, 
we restrict in Section [Hour description of the surface plas- 
mon to its two lowest eigenstates and discuss the resulting 
dynamics. The regime of validity of the two-level approxi- 
mation will be discussed in the appendix. In Section[5j the 
model parameters for typical experimental situations are 
estimated. Finally, we present in Section [S] our prediction 
of sidebands in the absorption spectrum of a weak probe 
beam impinging on a nanoparticle subject to a pump laser. 



2 Model of a laser-driven nanoparticle 

The electron gas formed by the N valence electrons of the 
metallic nanoparticle is described by the Hamiltonian 
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where is the position of the i th electron. The single- 
particle confining potential U represents the interaction 
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of the electrons with the ionic background which within 
the jcllium model is approximated by a uniformly charged 
sphere of radius a [2] . As a function of the radial coordi- 
nate, this central potential is harmonic for |r| < a and 
Coulomb-like for |r| > a. 

For our purposes it is useful to express the Hamiltonian 
(fT|) in terms of the centre-of-mass and relative coordinates 
as [11,13] 

H = H cm + H ie ] + H c . (2) 

For small centre-of-mass displacements, H cm is well de- 
scribed in terms of the usual ladder operators and b by 
the harmonic oscillator Hamiltonian 



H cm = M> M tfb 



with a frequency 



U)M\ 1 - 



N 
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Here, ujm = (47rn e e 2 /3m e ) 1 / 2 is the classical Mie frequen- 
cy which depends on the charge e, the mass m e and the 
bulk density n c of the electrons. The spill-out effect [2], 
which accounts for the fact that a fraction N out /N is found 
outside the nanoparticle, reduces the electron density and 
thus leads to a redshift of the plasmon frequency from the 
Mie frequency to wm- Except for extremely small clusters, 
anharmonicities yield negligible corrections to ((3]) and are 
thus disregarded here [15]. 

Within the mean-field approximation, the Hamiltonian 
for the relative coordinates reads 



(5) 



where e a are the eigenenergies in the effective mean-field 
potential V and (c a ) are the creation (annihilation) op- 
erators associated with the corresponding one-body eigen- 
states | a). Numerical calculations within the local density 
approximation show that, for the purpose of analytical cal- 
culations, V can in general be approximated by a step-like 
potential at the surface of the nanoparticle [12, 13]. 

The linearised coupling Hamiltonian between the cen- 
tre-of-mass and relative-coordinate systems can be written 
as 

H c = A (6+ + b) d aP clcp , (6) 

a/3 

where the coupling strength A — (hm c uj'l l /2N) 1 / 2 . The 
matrix element d a p between two eigenstates of the un- 
perturbed mean-field problem appearing in (JSJ) is for a 
spherical hard-wall potential simply given by the dipole 
matrix element (a|z|/3) [12]. Here, we have assumed that 
the surface plasmon oscillates in the z direction. The sum 
of ^ has to be understood with a low-energy cutoff since 
the electron-hole excitations of low energy do not couple 
to the surface plasmon [13,16]. This cutoff energy has been 
estimated in [17] to be of the order of 3/5 of the Fermi 
energy for the case of Na nanoparticles. 



The decomposition is reminiscent of the well-stu- 
died case where the degree of freedom of interest (the sur- 
face plasmon in our case) is weakly coupled to a reservoir 
with many degrees of freedom (the relative coordinates). 
Interestingly, it can be shown that it is sufficient to have 
about N = 20 conduction electrons to define a proper 
environment for the surface plasmon [17]. Since the elec- 
tronic system is coupled to phonons, a time dependence 
of the electronic temperature T may have to be taken 
into account. However, the physical parameters describ- 
ing the dynamics of the surface plasmon are only weakly 
temperature dependent [13], and thus we can neglect the 
implicit time dependence of those parameters when T is 
much smaller than the Fermi temperature of the system. 

If the nanoparticle is subject to a laser excitation, we 
have to add the coupling H-p (t) between the electrons and 
the laser field to the Hamiltonian ^ . Excitation frequen- 
cies wl close to the plasmon resonance are located in the 
visible range, and therefore the electromagnetic field has 
a wavelength much larger than the size of the nanopar- 
ticle. This external electrical field can be considered as 
spatially homogeneous, E(t) = E cos {u}\,t)e Z) and thus it 
only couples to the electronic centre of mass. In the dipolar 
approximation, the interaction between the driving field 
and the centre-of-mass system reads 



H F (t) = (6 f + b)HS7 R cos (uj h t) . 
The Rabi frequency i?R is given by 



hH R = eE Nl , 



where 



2Nm c CjM 



(J) 



(8) 



(9) 



is the characteristic length associated with the harmonic 
oscillator 

Different physical regimes can be attained according to 
the relative values of the typical energies Swm, aA, hf2 Rl 
and the Fermi energy of the metal. 



3 Centre-of-mass density matrix 

In the first part of this section we consider the free evo- 
lution of the centre of mass without the external driving 
field. The driving term ((TJ will be added in the second 
part. 



3.1 Free evolution of the centre of mass 

The evolution of the whole system in the absence of the 
electromagnetic field (Hp — 0) can be expressed in terms 
of the total density matrix W(t). Assuming that the cou- 
pling ([S]) is a small perturbation, we write the time evolu- 
tion of the density matrix as 



W(t) 



H c (t),W(t) 



(10) 
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where the tilde denotes operators in the interaction picture 
with respect to the uncoupled Hamiltonian H cm + H le i. 

As we are only interested in the dynamics of the centre 
of mass, we trace out the relative coordinates to obtain 
the reduced density matrix p = Tr^W^ which obeys the 
equation of motion 



Pit) 



n 



H c {t),W(0) 



1 



ds Tr rc 



H c (t), H c (s),W(s) 



(11) 



In dissipative quantum dynamics the environment is 
thought of as being infinite and remaining in thermal equi- 
librium at all times, i.e., its intrinsic properties like the 
temperature do not change as a consequence of its cou- 
pling to the system [18]. In order to satisfy this condition, 
we have to assume that the internal equilibration of the 
relative degrees of freedom after an excitation occurs suf- 
ficiently fast. This last assumption is not equivalent to 
the thermalisation of the whole electron gas, which can 
only occur when the plasmon has decayed. The internal 
thermalisation hypothesis can be justified a posteriori by 
checking that the timescale of the temperature change of 
the relative degrees of freedom is larger than all the other 
relevant timescales. 

Restricting ourselves to the lowest order in the corre- 
lations between system and environment in pip , the weak 
coupling limit that we assume to be valid allows to write 
(for details, see Chap. 4 in [14]) 



W{t) M p{t) <g> Ple} 



at any time t 0. Here, 



Prol 



(12) 



(13) 



represents the density matrix of the environment in the 
grand-canonical ensemble, i.e., we ignore for this purpose 
the finite nature of the electron gas. S is the grand-ca- 
nonical partition function at the inverse temperature (3 = 
(fceT) -1 and p is the chemical potential. 

As a consequence of the structure of the coupling H c , 
the fact that d aa — due to the dipole selection rules [12], 
and the assumption of a factorising initial condition, the 
first term on the right-hand side of (fTTj) vanishes. With 
the weak-coupling limit, (jTTJ) then reads 



m * - 



ds Tr 



rel 



H c (t), H c (s),p(s) <g> /9 re l 



The trace over the electronic environment §5§ yields 



(14) 



P® = ¥J dT C(T) 

6 + (f) + b(t),p{t - r) fo(t -t) + b(t - r)) 



h.c. 

(15) 



The correlation function of the environment C(r) is de- 
fined as 

C{t) = A 2 J2 [1 - /Ml /M Kef e^ r , (16) 

a/3 

and contains the relevant information on the time evo- 
lution of the relative-coordinate degrees of freedom. In 
(fig]) , f(e) = [e^ s ~^ + l]^ 1 is the Fermi function and 

For not too small nanoparticles, the electronic environ- 
ment described by the Hamiltonian ((5|) contains a large 
number of degrees of freedom, and its spectrum is there- 
fore quasi-continuous. Thus, the exponentials contributing 
to the correlation function (fl6|) very efficiently suppress 
the sum for not too small r. Indeed, for zero temperature 
it is shown in [17] that the typical correlation time (r cor ) 
of the environment in terms of the surface plasmon re- 
laxation time T p i is given by (T cor )/r p i « (fcpa) -1 , where 

is the Fermi wavevector. For nanoparticles with a ra- 
dius larger than approximately 1 nm, one typically finds 
r p i to be of the order of or larger than 10(r cor ). The effect 
of finite temperatures on this estimate is expected to be 
weak [13]. 

For sufficiently large nanoparticles, we can therefore 
work within the Markov approximation [14], where the 
density matrix p in (|15p can be assumed to vary on much 
longer timescales than the decay time of the correlation 
function C(t). The integral in (fT5|) can then be carried out. 
After returning to the Schrddinger picture, and within the 
secular approximation [14], where highly oscillating terms 
are neglected, we finally obtain the master equation for 
the reduced density matrix of the centre-of-mass degree 
of freedom in the absence of the external driving field 

p(t)=-iu sp [tfb,p(t)] 

- \ [b f bp(t) + p{t)tfb - 2bp(t)tf] 

-^[btfp(t)+p(t)btf-2tfp(t)b] . (17) 

Here, we have defined 

7± = tt^ 2 53 1 1 - f( £ *)\ /M \d Q fi\ 2 S (wm ± U) a/3 ) , 

a/3 

(18) 

while the renormalised surface plasmon frequency reads 
^sp = £>m - S (19) 

with 



s = i A2v E t 1 - /( £ «)] I'M 2 ■ 1 20 • 

a/3 



where V denotes the Cauchy principal value. 

The expressions for 7± and S entering (|17[) have simple 
physical interpretations [13]: Using Fermi's golden rule, 
7+ and 7_ are related to the lifetime r„ = j^ 1 of the 
n th excited state of the harmonic oscillator describing the 
centre-of-mass system by 



7n = "-7- + ( n + 1)7 



+ 1 



(21) 
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while S is the frequency shift due to the interaction of 
the system with the electronic environment. As the sur- 
face plasmon corresponds to the first excited state of the 
harmonic oscillator, its lifetime is given by r p i = j^ 1 . We 
have evaluated these important and experimentally rele- 
vant quantities in references [12] and [13] as a function of 
the size and temperature of the nanoparticle. The shift 
S is positive, and thus we have a redshift of the surface 
plasmon frequency in addition to the redshift induced by 
the spill-out effect described by (j4]). This redshift of the 
surface plasmon frequency is due to the coupling of the 
centre-of-mass system to the electronic environment. It 
is analogous to the Lamb shift known in atomic physics 
[14,19]. 



3.2 Effect of the laser 



We now consider the role of an external driving field whose 
interaction with the centre-of-mass coordinate is described 
by the time-dependent Hamiltonian Hp (t) denned by ([7]) . 
In a first step, we neglect the coupling Hamiltonian H c be- 
tween the centre-of-mass and the electronic environment 
during the excitation process. Thus, we do not have any 
relaxation mechanism for the surface plasmon excitation, 
and its frequency (|19|) remains unrenormalised. As a con- 
sequence, we have 

p(t) = — kI>m [tfb, p{t)} - if2 R cos (tu L t) [6 f + 6, p(t)] (22) 

which describes the standard time evolution of a harmonic 
oscillator driven by an external monochromatic field. 

Assuming that the driving does not influence the dis- 
sipation of the surface plasmon excitation, we can add the 
contributions from H c , i.e., the dissipative part of the re- 
duced density matrix described by the master equation 
(fT7)l . and from H F (t) (see Eq. ([22"])) independently. This 
is justified provided that the memory time of the environ- 
ment (r cor ) is much smaller than fi^ 1 which is determined 
by the coupling to the driving field [14]. The matrix repre- 
sentation of the master equation in the harmonic oscillator 
basis is then that of a driven damped harmonic oscillator, 



Pn 



iu S p(n - m)pr, 



( 12 + TO r — ; \ 

7 I ^ Pnm ~ V\ n + 1 ){ m + l)P»+l,m+l I 



iJ? R cos (w L i) (Vrip„-i,m + Vn+ lp„+i,, 



mp n , m -i - \/mTT/)„, m+ i) 



(23) 



In this equation, we have neglected 7+ as compared to 7_ 
and set 7 = 7- which corresponds to the Landau damping 
lincwidth. Indeed, it is easy to show from (JTSJ) that the rate 
7+ is related to 7_ through the detailed-balance relation 
7 + = e ~' 3;iWM 7_. We have 7+ <C 7_ for temperatures up 
to a few thousand degrees since Huim is of the order of 
several eV. 



4 Two-level system approach 

In order to obtain from (|23|) the time evolution of the 
centre-of-mass degree of freedom under the influence of 
the external driving field and the coupling to the relative- 
coordinate system, it is useful to introduce an appropriate 
simplification. The centre-of-mass system, which has been 
modelled as a harmonic oscillator described by the Hamil- 
tonian ([3]), can be truncated to a two-level system. The ap- 
plicability of such an approximation is discussed in the ap- 
pendix but can be motivated as follows: Except for a very 
strong driving field, the harmonic oscillator states above 
the first excited state are not significantly populated. Fur- 
thermore, the detuning between the frequency of the laser 
and the resonance frequency of the system plays in favour 
of the two-level description. Moreover, there exist addi- 
tional damping mechanisms which are not included in our 
model which tend to depopulate the higher states, like the 
ionisation of an electron via the double plasmon state [12] 
or the radiation damping [2] . In the case of pump-probe 
experiments on noble-metal nanoparticles embedded in a 
dielectric medium, interactions with the surrounding ma- 
trix provide further decay channels, e.g. via coupling to 
phonons, localised states, surface states, etc. As the second 
collective level, the so-called double plasmon, has a width 
which is significantly larger than the one of the simple sur- 
face plasmon, it is justified to neglect all excited levels but 
the first one. In addition, the double-plasmon state has not 
been clearly identified experimentally, even though indi- 
rect observations of such a state have been reported in 
experiments on charged sodium clusters in vacuum [20]. 

Writing the master equation (|23p for the two collective 
states |0) and |1), introducing the new variables 

Q — o piwL (n-m)t (2A) 

and keeping only the terms which significantly contribute 
close to the resonance cjl ~ ^sp (rotating wave approxi- 
mation [14]), one obtains 



pn 



Mr 

!— (Pio 



P01) - 7Pn > 



.Or 



P01 = -iS h poi - i— (pn 



Poo J 



7 



7P01 



(25a) 
(25b) 



Here, <5l = wl — <^sp is the detuning between the laser 
and the resonance frequency, and we have the conditions 
Poo + P11 = 1 and p\q = poi- It is often useful to define a 
scaled detuning A — <5l/7 and the saturation parameter 



7 



(26) 



which is a measure of the ratio between the field intensity 
and the damping strength. As a function of these param- 
eters, whose experimentally relevant values are discussed 
in the next section, the stationary solutions of the Bloch 
equations (I25|) are [14] 

(27) 



Pn — 



P01W 



2(1 



i + 4A 2 ) ' 
2A + i 
1 + s + AA 1 



(I) 



1/2 



(28) 
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The occupation probability pf* of the first excited state 
of the centre-of-mass system increases with increasing sat- 
uration parameter s and decreases with increasing detun- 
ing A between the frequency of the laser field ojl and 
the resonance frequency Lo sp . The coherence poi deter- 
mines the mean centre-of-mass coordinate according to 
(Z) = 2£Re( / ooi) where i has been defined in ((9|. The am- 
plitude of the oscillation of (Z) increases with increasing 
s up to s = 1 + 4Z\ 2 , while it decreases for larger s and 
approaches zero in the limit of infinite s. 

For zero detuning, <5l = 0, (f2"5")) can be solved analyti- 
cally. Assuming that the system is initially in its ground 
state, p(0) = |0)(0|, the full time dependence of the den- 
sity matrix is determined by [21] 



Pu(t) 



2(1 + s 
3 



1 



-37*/4 



sinh ( - 8s 

for the population of the excited state, and 

1/2 



cosh ( ^Vl - 8s 



(29) 



«*® = 1+7(2) 
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-37t/4 



cosl) ( ^j-Vl — 8s 



1 - 2s 



- ; sinh —VI - 8s 

^/T^8i V 4 



It 



(30) 



for the coherence. 

Figure [TJl shows the population (|2"9"|) of the excited 
state pn as a function of time for various values of the 
saturation parameter and without detuning (#l = 0). For 
s ^ 1/8, pn monotonically increases as a function of time 
to reach the stationary value given in (|27|) . For s > 1/8, 
the centre-of-mass system exhibits damped Rabi oscilla- 
tions and reaches the steady state on the timescale ~ 7 _1 . 
If the laser field is turned off after a certain time r, the 
population of the excited state relaxes to zero according 
to 

Pn (t) = Pn (T)e-^- T \ (31) 
while the coherence of the system decays as 

-( 7 /2-mj sp )(t-T) 



Poi(i) = Poi(r)e~ 



(32) 



One can see in Figure [T] that even for a strong saturation 
parameter, the steady state regime is reached for times 

In Figure [TJd, we show the mean centre-of-mass posi- 
tion (Z) as a function of time. For s <C 1, (Z) oscillates 
regularly, while its amplitude increases monotonically as 
a function of time. In the case of strong saturation param- 
eter s 3> 1 , it is clear from (jST)]) that the dynamics of the 
centre of mass is determined by the beating between the 
two oscillating contributions at the frequencies uj sp and 
J?rc for times shorter than a few 7 -1 (see the solid line in 
Pig. ID)). 



5 Saturation parameter 

The results in the previous section show that the satu- 
ration parameter (|26p is decisive for the dynamics of the 
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Fig. 1. (a) Population of the excited state pw (|29[) as a func- 
tion of 7t for different values of the saturation parameter s and 
for <5l = 0. (b) Mean centre-of-mass coordinate (Z) for s — 0.1 
and 100. In the figure, w sp /7 = 30 and 5l = 0. 



surface plasmon excitation in the presence of an external 
driving field. In typical photoabsorption experiments [1], 
a weak laser field excites an ensemble of nanoparticles. For 
a laser power of a few mW, the resulting electric field is of 
the order of 10 3 Vm~ 1 , and thus the Rabi frequency (JSJ 
is no more than several fieV. The order of magnitude of 
the surface plasmon linewidth in metallic nanoparticles is 
100 meV, so that one finds s <C 1. In that case we are in 
the linear-response regime, i.e., the power absorbed by the 
nanoparticle is proportional to the intensity of the laser. 

On the contrary, in pump-probe experiments, the na- 
noparticles are excited by an ultrashort pump laser pulse 
of high intensity. As a typical example, we consider the 
experiments of [7] on silver nanoparticles with an average 
radius a — 3.25 nm embedded in a glass matrix. These 
nanoparticles present a broad absorption spectrum around 
the resonance frequency aj sp = 2.85 eV. According to time- 
dependent local density approximation calculations [12], 
the width of this resonance is approximately 7 = 50 meV. 
The estimation of the Rabi frequency |(5J) requires the 
knowledge of the electrical field Eq of the laser. This quan- 
tity can be inferred from the energy u — e^E^cr /2 of a 
laser pulse divided by its cross section [7], where c is the 
speed of light, eo the vacuum permittivity, e r ~ 6 for glass, 
and r the duration of the pulse (typically 100 fs). With 
u = Urn" 2 , one finds E = 3.6 x'lO^uT 1 . The Rabi 
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frequency is therefore ft R ~ 0.38 eV, which results in a 
very large saturation parameter s ~ 100. Thus we are far 
beyond the linear-response regime. Notice that the dura- 
tion of the pump-laser pulse considered here corresponds 
to 7T ~ 8. Thus, according to Figure [TJ even for ultrashort 
laser pulses the centre of mass can reach its steady state 
before the end of the pulse. 



6 Absorption spectrum of a weak probe beam 

As it is shown in the appendix, the dynamics of the surface 
plasmon is well described by a two-level system driven by 
an external field and damped by its coupling to the elec- 
tronic environment. This is reminiscent of the well-known 
situation in quantum optics, where a transition between 
two levels of an atom is excited by a strong laser. The con- 
cept of dressed levels [14,22], arising as coherent superposi- 
tions of atom-photon states has proven to be very useful to 
describe the coupled system, in particular when the driv- 
ing is strong as compared to the damping. In the presence 
of a small detuning <5l between the pump laser field and 
the transition frequency of the two-level system, the ab- 
sorption spectrum of a second laser acting as a weak probe 
contains sidebands. Interestingly, in one of the sidebands 
the probe beam is damped by absorption processes, while 
stimulated emission processes amplify the probe beam in 
the other sideband. This phenomenon has been observed 
experimentally for different driven systems. Examples in- 
clude the case of an optical transition in sodium atoms [23] 
and a nuclear magnetic hyperfine transition of a nitrogen- 
vacancy centre in diamond [24] . 

In the present paper we focus on a surface plasmon ex- 
citation driven by a strong laser field. This is a completely 
analogous situation since the dynamics of the surface plas- 
mon can be described by an effective two-level system 
with a transition energy huj sp , coupled to a pump laser 
with frequency Without coupling = 0), the basis 
states are \0,N P + 1) and 1 1 , iV p ) representing the plas- 
mon in its ground and excited state with N p + 1 and N p 
photons in the pump laser field, respectively. The energy 
difference between those states is given by the detuning 
Ml = K(utL — u S p) of the pump laser with respect to the 
plasmon transition. The coupling by Hp (Eq. 0) leads 
to a mixing of the basis states, giving rise to the dressed 
plasmon levels with an energy splitting HQ — K^J fti + 8^, 
where J?r is the Rabi frequency fH} . The dressed plasmon 
levels \A(N P )) and \B(N P )) arise as superpositions of the 
uncoupled basis states |0, N p + 1) and [1, JV P ) referring to 
different photon numbers. In this doublet, we denote the 
upper and lower dressed level by \A(N p )} and \B(N P )), 
respectively. 

Four optical transitions (see thick arrows in Fig. [5]) are 
allowed between dressed states corresponding to neigh- 
bouring photon numbers because of the mixture of those 
states with the ones of the uncoupled basis. Since the sta- 
tionary populations of \A(N p )) and \B(N p )) are to a very 
good approximation independent of N p , there is no peak 
in the absorption spectrum of the weak probe beam at the 
frequency ojl (see the two thick dashed arrows in Fig. [2]). 



\l,N p ) 



|0,JV P + 1> 



|1,JV P -1> 



|0,iVp> 




B(N P )) 



A(N P - 1)) 



B(N P - 1)) 







Fig. 2. Scheme of the uncoupled (4?r = 0) and dressed states 
(J?r 7^ 0). The four allowed optical transitions are indicated 
by vertical thick arrows. See the text for details. 



For Si, > 0, the upper dressed levels \A(N P )) have a 
higher stationary population than the lower ones \B(N P )). 
Therefore, a population inversion appears in the transition 
\A(N P )) <-> \B(N P -1)). This leads to a sideband at the fre- 
quency wl + ft where the probe beam is amplified by stim- 
ulated emission. For the transition \B(N P )) <-> \A{N V — 1)) 
the lower state has the higher population such that the 
sideband at frequency u^—fi is of absorbing character. For 
negative detuning <5l < 0, when the pump laser frequency 
is below the bare transition energy, the populations of the 
dressed states are inverted. Now, the populations of the 
upper dressed levels \A(N p )) are lower than the popula- 
tions of the lower levels. In this case, the sideband at the 
frequency wl + ft is absorbing while the probe beam is 
expected to be amplified around the frequency wl — ft- 

The width of the sidebands is of the order of 7. In the 
so-called secular limit ft ^> 7, the sidebands are there- 
fore well resolved and should be observable experimen- 
tally. Within this limit, it is straightforward to calculate 
the stationary populations of the dressed states [22]. For 
the case of a resonant pump laser 8^ — 0, all the dressed 
states have the same populations. Therefore, to zero order 
in j /ft, the transitions do not influence the probe beam. 
A more detailed discussion beyond the parameter regimes 
considered here is given e.g. in references [23] and [25]. 

An experimental observation of the sidebands should 
be feasible using an intense pump laser with s ~ 100 and 
a pulse duration longer than 100 fs. Then, the frequency 
width of the pump pulse is below 0.01 eV, which allows to 
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adjust the detuning Sl with high precision. With the Rabi 
frequency Q-r — 0.38 eV and the lincwidth 7 — 50meV 
estimated in Section [5] for the experiments of reference [7], 
one should be able to observe well resolved sidebands. 



7 Conclusion 

Using standard methods of quantum optics, we have stud- 
ied the dynamics of the surface plasmon in metallic nano- 
particles driven by an electrical field. Decoherence and dis- 
sipation occur due to the coupling of the collective excita- 
tion to the electronic environment constituted by particle- 
hole pairs. Exploiting the Markovian character of the en- 
vironment [17], we have established the master equation 
for the surface plasmon density matrix. 

We have shown that for realistic experimental situ- 
ations, a model taking into account only the two lowest 
levels is sufficient to describe the dynamics of the collective 
degree of freedom. This allows to consider a large variety 
of experimentally relevant situations, from photoabsorp- 
tion to pump-probe setups. 

Motivated by the analogy to quantum optics, we have 
introduced collective states that are dressed by the pho- 
tons of the pump laser. We predict that for strong driv- 
ing, the dressed states will manifest themselves through 
the appearance of sidebands in the absorption spectrum 
of a weak probe laser. Depending on the detuning of the 
pump laser, these spectral features have absorbing or am- 
plifying character and their position relative to the pump 
frequency is determined by the pump intensity. Although 
these sidebands have so far not been observed in metal- 
lic nanoparticles, their detection should be feasible with 
modern experimental techniques, thereby offering access 
to dressed surface plasmon states. 
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phase of this work. 



A Three-level system 

In this appendix we discuss the validity of a two-level ap- 
proach for the description of the centre-of-mass system, 
whose Hamiltonian is given in ([3]). For this purpose, we 
consider a three-level system and study the conditions un- 
der which its description can be safely reduced to that of 
a two-level system. A phenomenological master equation 
description of the plasmon in terms of three levels has been 
discussed by Liau et al. [26]. In that case, however, three 
levels where required to describe two-photon processes in 
two-pulse second-order interferometry. 

According to (|2"T|) . the Landau damping linewidth for 
the second level of the harmonic oscillator at zero tem- 
perature is 27, where 7 is the width of the first excited 



state [12,27]. In addition, two more mechanisms contribute 
to the damping of the second level: (i) first order processes 
which lead to the decay into the first level (like the radi- 
ation damping [1,2] or the coupling to the surrounding 
matrix [27]); (ii) second-order processes which result in 
the direct decay into the ground state (like the ionisa- 
tion [12]). We denote the damping rates associated with 
these two additional channels by 71 and 72, respectively. 

Within the Lindblad theory [28], these channels are 
accounted for by adding 

(n\ J2 y {2L a p4 - L\L aP - ph\L a ) \m) . (33) 

a=l,2 

to the right-hand side of the master equation (|23|) . We 
choose the Lindblad operators as L\ = |1)(2| and L 2 = 
|0)(2|, which lead to transitions between the centre-of- 
mass states 1 2) and |1), and |2) and |0), respectively. The 
rotating-wave approximation then yields the following set 
of coupled differential equations 



Poo = 
P22 = 

P01 = 
P12 = 



Mr 



1— (pio — P01) + IPll + 72P22 ■ 
i£|iV2(/5i2 - P21) 



(27+ J2 th 

\ 0=1,2 / 



P22 



i^LPoi - i-^p (pn - Poo - V%po2 



:Poi 



V2jpi2 , 

\P02 



P02 = 



iShPl2 - i-^ 1 [p02 + V2 (p 2 2 - Pll) 

2 ( 37 + Y 7o J Pl2 , 

V o=l,2 / 

2WlPo2 - i-^ 1 (pi2 - Vlpoi 

\\ 27 + 7q ) P° 2 ■ 

\ 0=1,2 / 



(34) 



In order to obtain quantitative results, the two addi- 
tional damping constants in the following are assumed to 
be equal, i.e., 71 = 72 = 7add- Figure [3] presents the sta- 
tionary solution of (f34|) for the population of the second 
excited state p2% as a function of 7 a dd for different values of 
the detuning A = <5l/7 and a large saturation parameter 
s = 100. It can be seen that for A > 5, which corresponds 
to the typical detuning in the pump-probe experiments of 
reference [7], the probability to find the second collective 
state occupied is less than 10% and is almost constant as 
a function 7 a dd- In the case A < 5, pf 2 decreases with 
7add but may be non-negligible. In this case, the assump- 
tion of a two-level system to describe the surface plasmon 
dynamics might be questionable. 

For a rather small saturation parameter, i.e., for a weak 
external laser field, the second collective state can be ne- 
glected for the description of the surface plasmon dynam- 
ics. In the case of a large saturation parameter, as it is for 
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Fig. 3. Stationary solution of (|34p for the population of the 
second excited state as a function of the additional damping 
constant 7 a( jd and for various values of the detuning A. In the 
figure, the saturation parameter is s = 100. 



example the case in pump-probe experiments, this second 
collective state can be neglected for not too small detun- 
ing between the frequency of the pump laser and the res- 
onance frequency. 



18. U. Weiss, Quantum Dissipative Systems (World Scientific, 
Singapore, 1993) 

19. W.E. Lamb and R.C. Retherford, Phys. Rev. 72, 241 
(1947); W.E. Lamb, Rep. Prog. Phys. 14, 19 (1951) 

20. R. Schlipper, R. Kusche, B. von Issendorff, and H. Haber- 
land, Phys. Rev. Lett. 80, 1194 (1998); Appl. Phys. A: Mater. 
Sci. Process. 72, 255 (2001) 

21. D.F. Walls and G.J. Milburn, Quantum Optics (Springer- 
Verlag, Berlin, 1994) 

22. C. Cohen- Tannoudji and S. Reynaud, J. Phys. B: Atom. 
Molec. Phys. 10, 345 (1977); C. Cohen-Tannoudji and S. 
Reynaud, in: Proc. Int. Conf. on Multiphoton Processes, 
ed. by J.H. Eberly and P. Lambropoulos (Wiley, New York, 
1978) 

23. F.Y. Wu, S. Ezekiel, M. Ducloy, and B.R. Mollow, Phys. 
Rev. Lett. 38, 1077 (1977) 

24. C. Wei and N.B. Manson, Phys. Rev. A 49, 4751 (1994) 

25. B.R. Mollow, Phys. Rev. A 5, 2217 (1972) 

26. Y.-H. Liau, A.N. Unterreiner, Q. Chang, and N.F. Scherer, 
J. Phys. Chem. B, 105, 2135 (2001) 

27. R.A. Molina, D. Weinmann, and R.A. Jalabert, Phys. Rev. 
B 65, 155427 (2002); Eur. Phys. J. D 24, 127 (2003) 

28. G. Lindblad, Commun. Math. Phys. 48, 119 (1976); V. 
Gorini, A. Kossakowski, and E.C.G. Sudarshan, J. Math. 
Phys. 17, 821 (1976) 



References 

1. W.A. de Heer, Rev. Mod. Phys. 65, 611 (1993) 

2. M. Brack, Rev. Mod. Phys. 65, 677 (1993) 

3. T. Tokizaki, A. Nakamura, S. Kaneko, K. Uchida, S. Omi, 
H. Tanji, and Y. Asahara, Appl. Phys. Lett. 65, 941 (1994) 

4. J.-Y. Bigot, J.-C. Merle, O. Cregut, and A. Daunois, Phys. 
Rev. Lett. 75, 4702 (1995) 

5. M. Perner, P. Bost, U. Lemmer, G. von Plessen, J. Feld- 
mann, U. Becker, M. Mennig, M. Schmitt, and H. Schmidt, 
Phys. Rev. Lett. 78, 2192 (1997) 

6. V. Halte , J. Guille, J.-C. Merle, I. Perakis, and J.-Y. Bigot, 
Phys. Rev. B 60, 11738 (1999) 

7. J.-Y. Bigot, V. Halte, J.-C. Merle, and A. Daunois, Chem. 
Phys. 251, 181 (2000) 

8. N. Del Fatti, F. Vallee, C. Flytzanis, Y. Hamanaka, and A. 
Nakamura, Chem. Phys. 251, 215 (2000) 

9. C. Voisin, N. Del Fatti, D. Christofilos, and F. Vallee, J. 
Phys. Chem. B 105, 2264 (2001) 

10. G. Weick, D. Weinmann, G.-L. Ingold, and R.A. Jalabert, 
Europhys. Lett. 78, 27002 (2007) 

11. L.G. Gerchikov, C. Guet, and A.N. Ipatov, Phys. Rev. A 
66, 053202 (2002) 

12. G. Weick, R.A. Molina, D. Weinmann, and R.A. Jalabert, 
Phys. Rev. B 72, 115410 (2005) 

13. G. Weick, G.-L. Ingold, R.A. Jalabert, and D. Weinmann, 
Phys. Rev. B 74, 165421 (2006) 

14. C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-photon interactions: Basic processes and applications 
(Wiley- VCH, New York, 1992) 

15. K. Hagino, Phys. Rev. B 60, R2197 (1999) 

16. C. Yannouleas and R.A. Broglia, Ann. Phys. (NY.) 217, 
105 (1992) 

17. C. Seoanez, G. Weick, R.A. Jalabert, and D. Weinmann, 
|arXiv:cond-mat /0703720 



